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Nuclear spins are promising candidates for quantum information processing because their good 
isolation from the environment precludes the rapid loss of quantum coherence. Many strategies have 
been developed to further extend their decoherence times. Some of them make use of decoupling 
techniques based on the Carr-Purcell and Carr-Purcell-Meiboom-Gill pulse sequences. In many 
cases, when applied to inhomogeneous samples, they yield a magnetization decay much slower than 
the Hahn echo. However, we have proved that these decays cannot be associated with longer 
decoherence times as coherences remain frozen. They result from coherences recovered after their 
storage as local polarization and thus they can be used as memories. We show here how this 
freezing of the coherent state, which can subsequently be recovered after times longer than the 
natural decoherence time of the system, can be generated in a controlled way with the use of field 
gradients. A similar behaviour of homogeneous samples in inhomogeneous fields are demonstrated. 
It is emphasized that the effects of inhomogeneities in solid state NMR, independently of their origin, 
should not be disregarded as they play a crucial role in multipulse sequences. 
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I. INTRODUCTION 

The initial proposals for quantum information 
processing^ triggered a worldwide effort to develop simple 
specific realizations* 2 * 3 -. In the last decade, much effort 
has been devoted to understand and control decoherence^ 
as it has become the main drawback for quantum infor- 
mation processing^—. 

Within solid state Nuclear Magnetic Resonance the in- 
terest focuses in those systems which have been proposed 
as possible candidates for quantum registers— ~—, quan- 
tum channels— ~—, simulators of quantum dynamics 17 i 18 
or probes for specific quantum gates^. 

More, recently, a number of interesting works seek to 
control the nuclear and electron spin entanglement in sys- 
tems where the nuclear spins can be used as quantum spin 
memories such as N@Cecr^, or 31 p^£ or nitrogen vacan- 
cies (NV centers) in diamond^— . Some of the strate- 
gies implemented with NMR and EPR, such as decou- 
pling techniques 2 ^ - — ., have a wide interest as they can 
be adapted to many other two level systems like trapped 
ions, quantum dots 2 ^ or Josephson junctions^— . The 
renovated interest in pulse sequences that help to fight 
decoherence has revived the old sequences such as the 
Hahn Echo2i, the Carr-Purcell (CP) & and the Carr- 
Purcell-Meiboom-Gill (CPMG)^ and it has stimulated 
the development of new one o 37 ' 38 . The simplest dynami- 
cal decoupling (DD) sequences consist on the application 
of 7r pulses on the system to revert the decay due to 
the system-environment interactions. Many works have 
been done in order to compare the performance of the 
CPMG-like sequences and de Uhrig (UDD) sequences^ 8 - 
when the system is coupled to a spin bat h 27 ' 39 . Yang, 
Wang and Liu have made a comprehensive theoretical 
and experimental description of the DD techniques and 
its variations^. 



In this work we will make a critical review of recent 
attempts to extend the decoherence time of spins in solid 
state NMR or EPR. We show that the role of line in- 
homogencity should not be underestimated to observe 
coherences reliably. 

We summarize our own findings on anomalous long 
lived echo signal o 12 ' 41 " — . In particular we emphasize the 
use of a simple stimulated echo sequence to discern if one 
is in the presence of true long decoherence times or if the 
coherences are being saved as polarization and then re- 
covered. This understanding allows us to develop a new 
strategy for their use that is reported here. It makes use 
of magnetic field gradients that enable to store part of 
the phase information as polarization. We demonstrate 
that inhomogeneous samples in presence of homogeneous 
fields behave in the same manner than homogeneous sam- 
ples in presence of inhomogeneous fields. 

Additionally, the correct phase cycle to use in multi- 
pulse sequences in order to measure the true T2 in in- 
homogeneous samples is shown. This sequence, whose 
result has been empirically known 4 ^— , should be taken 
as a standard to obtain a clean T 2 measurement. Here, 
we explain how the stimulated echoes are canceled out 
eliminating their main effects: the long tails. 



II. BASIC NMR CONCEPTS 

A typical NMR experiment requires at least two mag- 
netic fields: a static field Bo, conventionally in the z di- 
rection, and a perpendicular radiofrcquency (RF) field 
Bi applied to induce transitions. The time t p during 
which the RF is on is defined as the pulse time and the 
tilting angle is P = 7Bit p , where 7 is the gyromagnetic 
factor. The eigenstates of angular momentum for a spin 
1/2 along the static field direction are |t) and Once 



2 



the thermal equilibrium is reached the population of state 
|t) is larger than that of \\) indicating a net polarization 
of spins along the static field direction. This thermal 
polarization will determine the maximum magnetization 
that can be created in the system. Any other situation 
in which the macroscopic magnetization points in an ar- 
bitrary direction will eventually have a z-magnctization 
which reflects a non equilibrium polarization. A particu- 
lar but very relevant case is the application of a ^ or sat- 
uration pulse after which the magnetization lies in the xy 
plane. The population of the states is equilibrated and 
converted in a superposition state (a coherence) where 
the eigenstates |t) and \\) have the same probability and 
a well defined phase relation. The essential difference 
between polarization and coherences is that, in a very 
good approximation, polarization does not evolve under 
interaction with the static field while the transverse mag- 
netization (or coherence) does. 

Thus, after a ^ pulse is applied, transverse magneti- 
zations associated with different volume elements (called 
isochromats) evolve under different Hamiltonians due to 
small inherent inhomogeneities and spin-spin interac- 
tions. This evolution causes a dephasing among the spins 
in the sample and a free induction decay (FID) with char- 
acteristic time T2 is observed. The dephasing caused by 
inhomogeneities is reversible, so that refocusing is possi- 
ble by applying appropriate pulse sequences such as the 
Hahn echo^i [HE: (^ — r — n — t— echo)]. During the 
evolution time r a dephasing occurs and then a 7r pulse 
is applied. The net effect of this pulse is to change the 
sign of all the Hamiltonians linear in spin operators, thus 
reversing the evolution of the isochromats which refocus 
after the same time r. As a result a magnetization echo is 
observed, i.e. the coherent state created after the initial 
■| pulse is partially recovered after a time 2r. The ampli- 
tude of the echo obtained after a single refocusing pulse 
is attenuated by homonuclear interactions represented by 
non linear terms in the system Hamiltonian. In solids the 
most important contribution will be the spin-spin dipolar 
interaction. If the separation time r between the pulses 
is increased a decay in the magnetization is observed with 
a characteristic time T2. This quantity is known as the 
spin-spin decay time, transverse relaxation or dccoher- 
ence time. Actually, T2 represents the decoherence of a 
single spin caused by its interactions with other spins af- 
ter canceling out the interactions linear in spins. Indeed, 
the spin-spin interaction can also be reversed by suitable 
sequences. They extend the decoherence time to a new 
timescale T3 > T2 characterizing the ability to control 
the many-body dy namic a 1 7 ' 48 ~ ^£. Since, these sequences 
revert many-body interactions in recent years the quan- 
tum information community refers to them as Loschmidt 
Echoe s 17 ' 51 " —. In any case, the decoherence time T2 is 
the characteristic time that should be very well known 
before any attempt of quantum processing. It imposes 
a first natural limit to quantum control, beyond which 
more sophisticated experiments that restore phase infor- 
mation (e.g. decoupling and Loschmidt Echo time re- 



versal strategies) are neede d 17 ' 48 " —. Since 1950, several 
sequences have been developed in NMR to obtain the T2 
decay time. As the Hahn echo experiment is very time 
consuming, multipulse sequences that enable the acqui- 
sition of the full decay with a single train of pulses such 
as the Carr-Purcell^ 5 - CP : (?) s — [r — ir x — r — echo} n 
and the Carr-Purcell-Meiboom-Gill^ 6 . CPMG : (%) - 



echo 



are the most commonly used. 



III. LONG LIVED SIGNALS 

The- interest in nuclear spins for QIP, in particular sil- 
icon^ 5 -, focused on systems where nuclei couple to elec- 
tron spins which are good candidates for hybrid tech- 
nologies. This motivated many NMR works where typi- 
cal experiments to measure T2 were reporte d 12 ' 13 ' 56 . In 
these works after applying the CPMG like sequences it 
was found that it is possible to detect the spin- 1/2 NMR 
signal of 29 Si up to times much longer (more than two 
orders of magnitude) than the characteristic decay time 
observed with the Hahn echo sequence (« 5.6 ms). This 
striking and promising finding came as a surprise and 
triggered the realization of many experiments to under- 
stand its origin. From the beginning it was verified that 
there is no dependence on the amount of donors or accep- 
tors in the 29 Si polycrystalline sample even when these 
could change the linewidths and spin-lattice relaxation 
time by one order of magnitude. Experiments were re- 
peated on 29 Si single crystals and Pyrex, as well as 13 C 
in Ceo yielding the same behavior. 

Since the measurements were performed in natural 
abundance 29 Si (4.67%) and 13 C (1.1%), it was clear that 
the spins are diluted yielding weakly coupled networks. 
Notice that the spin dipolar interactions decay with the 
cube of the distance and that the gyromagnetic factors of 
29 Si and 13 C are approximately one fifth and one fourth of 
the X H rendering 29 Si- 29 Si and 13 C- 13 C couplings which 
arc 1/25 and 1/16 of the 1 H- 1 H respectively. This sug- 
gests that the disorder of the chemical shift energies (or 
site energies) in such weakly coupled networks could play 
a relevant role in the observed phenomena, which have 
not been observed in the popular proton samples (ie. in 
solid state NMR of : H the usual observation is that the 
results from a Hahn echo coincides with those of the CP 
or CPMG sequences). With this in mind we started a 
series of experiments and calculations in a natural abun- 
dance polycrystalline sample of Cgo- The first experiment 
was a traditional Hahn echo measurement as a function 
of the time r which yielded r Y^ E ~ 15 ms. It should be 
noted that the decay of the FID revealed a T2 ~ 2 ms, 
indicating a high degree of inhomogeneties or reversible 
single spin interactions. Then, we used a series of dif- 
ferent trains, derived from the CP and CPMG sequences 
These sequences were named respecting their origin (CP 
or CPMG for ir pulses with 0° or 90° phase shifts with 
respect to the first -| pulse, respectively) and considering 
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FIG . 1 : Comparison of the 13 C magnetizations acquired with 
different T2 sequences (r = 1 ms). The insert shows the decay 
observed with the Hahn echo and the short time regime of the 
CPMG2 sequence. 



FIG. 2: C Magnetization decay as a function of time mea- 
sured for different time windows with the CPMG2 sequence 
in polycrystalline C60 • It is evident that the envelope of the 
oscillatory decays are dominated by T2HB. 



the number of phases (1 or 2) for the ir pulses in each 
cycle as follows: 
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CPMG2 : [ - 
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t — (71") — r — echo 



(tt) — t — echo 
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Note that CP1 and CPMG1 arc the original CB^ and 
CPMG^i sequences respectively. 

We found that by applying the CPMG1 sequence or 
the CP2 sequence, a magnetization tail appears, i.e. the 
signal remains for times much longer than T^ E as ob- 
served in 29 Si,— . Figure Q] shows the normalized magne- 
tizations acquired by applying the sequences HE; CP1; 
CP2; CPMG1 and CPMG2 to polycrystalline C 60 . In 
the last four sequences r = 1ms and signal acquisition is 
performed at the top of the echo. In all the sequences the 
durations of the 7r and 7r/2 pulses were carefully set from 
nutation experiments to 7.0 fis and 3.4 fj,s, respectively. 

On the other hand, by applying the CP1 or the 
CPMG2 sequences the decay time is not longer than 
T^ E , but it is remarkable that at short times the magne- 
tization makes a zigzag (consequence of the presence of 
stimulated echoes) see inset in Fig|T] Besides, the mag- 
netization oscillates going through negative values before 
reaching its final asymptotic zero value but keeping the 
T^ E decay as envelope (see FigF2J). 

By plotting these oscillating decays as a function of the 
echo number a match of the frequencies can be observed, 
i.e. the maxima and minima overlap for the same num- 
ber of 7r pulses (echo number, see Fig|3]). This behavior 
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FIG. 3: C Magnetization decays as a function of the number 
of pulses (or echoes) for short r's in the CPMG2 sequence. 



reveals an accumulation of pulse errors which complete 
a cycle after approximately 80 pulses. The amplitudes 
decrease with increasing r following the behavior of the 
long tails. Note in FigQ]that for r = 1 ms,the oscillation 
almost disapears. 

Another unexpected observation was the presence of 
noticeable stimulated echoes [STE: (f ) 2 . — t — (§) — 
£1 — (?) ~t] when applying pulse sequences of the form, 



V2 



t\ — (t) — acq] with variable r 



and t\. Indeed, these sequences with ir pulses should 
not produce stimulated echoes. To remember the essen- 
tials of the argument we reproduce here our results af- 
ter the application of the STE sequence with the phases 
of the CPMG1 and CPMG2, called CPMG1 STB and 
CPMG2 ST£; respectively in Cgo setting r = 1 ms and 
t\ = 15 ms (see FigfJJ. It can be observed that the stim- 
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FIG. 4: Stimulated echoes observed with the CPMG1 STE and 
CPMG2 STB sequences in a Cgo polycrystalline sample. The 
parameters are r = 1 ms and t\ — 15 ms, giving rise to a 
stimulated and a normal echo at 1 ms and 14 ms after the 
last pulse, respectively. Notice the inversion of the phase in 
the stimulated echo. 



ulated echoes appear 1 ms after the last pulse (remember 
that they should not be formed with perfect it pulses) 
and they show a phase that coincides or it is opposite 
to that of the Hahn echo, that appears at t\ — r = 14 
ms, for the CPMG1 ST£; and CPMG2 STB respectively. 
A phcnomcnological explanation for the long magnetiza- 
tion tails based on the constructive or destructive inter- 
ferences between the stimulated and normal echoes was 
reported by us in a previous paper— and later on verified 
through hole burning experiments^ 4 -. Thus, although dis- 
sapointing for the experimentalists of the quantum infor- 
mation community, we concluded that the observed long 
tails in the magnetization decays are not a signature of 
long dccoherence times but of the recovery of coherences 
saved as polarization. These long pseudocoherent tails 
appear as a consequence of the formation of stimulated 
echoes. These echoes, which need at least three pulses 
to appear, were first observed by E. Hahn^i and have 
the particularity that can survive for times as long as 
the spin-lattice relaxation time Ti. They will be revised 
in the next section in order to try to save and recover 
coherences in a controlled way. 

This is the key experiment, and we suggest that this 
stimulated echo sequence with n pulses, [-| — t — n — 
It — 7r — t] should be used in the same conditions of 
the CP-like train sequences to be able to discard the 
stimulated echoes as the origin of the long magnetization 
tails. Although the presence of the stimulated echoes af- 
ter this sequence is not expected, they can be generated 
by the pulse angle distribution present in an inhomoge- 
neously broadened line. In all the samples mentioned 
above, the spin-spin decay time measured with a Hahn 
echo sequence (T 2 he)^ was about an order of magnitude 
longer than the FID characteristic time (Tj), evidencing 
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FIG. 5: Traditional stimulated echo sequence. The coherence 
states created during the hrst evolution period r, are stored as 
2-amplitude (polarization) during tiand then reconverted to 
coherences (transverse magnetization) after the third pulse. 



the line inhomogeneity. 

Magnetization tails have been observed in 29 Si, Cgo 
and Y20 ^ 13 ' 41 ' 42 , in adamantane under 1 H decoupling^ 
and indirectly in 31 P in EPR experiments^. 

The requirements for the formation of the long tails can 
be summarized as follows: (a) an rf field inhomogeneity 
or a naturally inhomogeneous line able to produce differ- 
ent tilting angles in different sites of the sample and (b) 
the absence of or a very slow spin diffusion (noneffective 
flip-flop interactions). Under these conditions the differ- 
ent tilting angles do not average to the nominal tt pulses 
preset in the sequences leaving a perpendicular compo- 
nent that behaves as in the regular stimulated echo se- 



quence 



2t 



acq 



schematized in FigJSJ Still, one question that remained 
open was if the relation between the spin-site energy dif- 
ference and the dipolar couplings (the main source of 
flip-flop interactions in solid systems) should be verified 
in a local (microscopic) basis or if a smoothly varying 
inhomogeneity that only satisfy the relation at longer 
distances could produce the effect. To investigate this 
issue we resorted to the use of gradients in homogeneous 
samples. 

We artificially create the necessary conditions to ob- 
serve the long tails after a CPMG measurement, even for 
homogeneously broadened samples. The method uses the 
fact that homogeneous lines, evolving under the influence 
of a field gradient, are inhomogeneously broadened giving 
rise to stimulated echoes with the consequent long tails. 
This method provides a controlled way to move from ho- 
mogeneous lines (for which the decay observed after a 
Hahn echo and CPMG like sequences are in complete 
agreement) to inhomogeneous ones. We reach inhomo- 
geneous limits analogous to 29 Si or Ceo cases, where we 
show that careful interpretation of the experiments must 
be done before assigning "the true" spin-spin decay time. 
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IV. BEHAVIOUR OF THE MAGNETIZATION 
DURING A SIMPLE STIMULATED ECHO 
SEQUENCE 

The stimulated echo contribution to the CPMG se- 
quence has been widely studied in presence of inhomo- 
geneous Bo and Bi fields. It has been demonstrated 
that due to the inaccuracy of the pulses and the off res- 
onance effects, the CPMG echoes have contributions not 
only from the Hahn echo but also from signals similar to 
stimulated echoes arising from different coherence path- 
ways^—. This is in complete agreement with our find- 
ings. 

Here, we review in detail the magnetization evolution 
during a stimulated echo to make evident how it can be 
potentially exploited to create memories. 

The cigenstates of angular momentum for a spin f/2 
along the static field direction are |t) and |4<). Once the 
thermal equilibrium is reached the population of state 
|t) is larger than that of |t) indicating a net polarization 
of spins along the static field direction. This thermal 
polarization will determine the maximum magnetization 
that can be created in the system. Any other situation 
in which the macroscopic magnetization points in an ar- 
bitrary direction will eventually have a z-magnetization 
which reflects a non equilibrium polarization. A particu- 
lar but very relevant case is the application of a ^ or sat- 
uration pulse after which the magnetization lies in the xy 
plane. The population of the states is equilibrated and 
converted in a superposition state (a coherence) where 
the eigenstates |t) and \l) have the same probability with 
well defined phase relation. The essential difference be- 
tween polarization and coherences is that, in a very good 
approximation, polarization does not evolve under inter- 
action with the static field while the transverse magneti- 
zation (or coherence) does. 

It is important to correctly understand the processes 
behind the stimulated echo in order to separate its contri- 
bution from the Hahn echo contribution. To clarify the 
details, let us go back to the simplest stimulated echo 
sequence: 
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In Fig. [5] the stimulated echo formation is schematized 
using a vectorial spin model. The model provides a picto- 
rial view where it can be observed that the magnetization 
is indeed stored as polarization and then reconverted to 
transverse magnetization. The information on the co- 
herence state (characterized by a well defined phase rela- 
tionship between the transverse spins) created in the first 
evolution period, t, is conserved during t\ and refocused 
as an echo during the third period at t = r. 

Let us consider a Hamiltonian linear in the spin oper- 
ators with no interactions between different spins, 



N 

H = Y,Wi- (3) 

3=1 

If the delta pulse approximation is taken and we con- 
sider the system in thermal equilibrium before the first 
pulse, the evolution of the density matrix can be written 
as: 

N 

p(r-) = - £ [cos(^-r)/J - sm(^r)/?] (4) 
j=i 

The expression above represents a coherent state 
(transverse magnetization), which immediately after the 
second pulse becomes, 

N 

P^ + ) = E [cos(^r)7J + sm(^r)/;] (5) 
3=1 

After the second evolution period (ti), the new state 
of the system can be expressed as: 

N 

p(r + ti) = — {cos(5wj-t)[/J cos(Sujjti) — I? sin(5ujjti)] 
+ sm^r)/?} (6) 

During the time ti, there is a component of the initial 
magnetization stored as non equilibrium polarization in 
the Ij. This component does not evolve during ti un- 
der the chemical shift Hamiltonian in Eq. ©. The last 
term within the sum in Eq. ([5]) defines a polarization grat- 
ing which lodges in its amplitude the information of the 
coherent state, i.e. the phase, of Eq. |QJ. The useful 
interpretation is that the coherent state created during 
the first evolution period remains frozen as a local po- 
larization amplitude during a long period t\. It should 
be noted that during the second period the polarization 
in Eq. ^ is not in thermal equilibrium (there is not a 
Boltzmann distribution because of the 5ui dependence). 
The spin-lattice relaxation, characterized by Ti, will lead 
to a decay of the polarization grating during t\. 

The third (§) pulse will convert the polarization of 
Eq. ([6]) in transverse magnetization as 

N 

p{r + tf) = ~y] {cos(<5wjt)[/J cos(5ujjti) + I] sin(£wjti)] 
3=1 

+ sin(^-r)/;} (7) 
which will evolve under the linear Hamiltonian influence 



6 



as: 

N 

p(t + r + ti) = — cos(8uijT) {cos(Sujjti)[iy cos(5uijt) 

3=1 

-If sin(5ujjt)] + I- sin(8wjti)} 

+ sm(SwjT) (ij cos(Sujjt) + Ij sm(5u)jt)) (8) 

The last two terms in Eq. <j8j) do not depend on t\ and 
when t = t the phase differences acquired between the 
spins in the first evolution period are refocused forming 
the stimulated echo. As can be seen from Eq.©, if ti=0 
the sequence is reduced to a simple (%) x — t — (w) — r— 
Hahn echo sequence and all the initial polarization is 
again at the -y axis. On the other hand, if t\ 3> T2HE, 
but still shorter than the spin-lattice relaxation time, 
only the term that do not depend on t\ survives and the 
coherences created during the first evolution period are 
refocused at t = r. Notice that if 2r ps T2HE the dipo- 
lar interaction can not be neglected and the assumptions 
performed in our calculations, where only the chemical 
shifts were considered, are not valid anymore. 

It is straightforward to see that if the second and third 
pulses in Fig. [5] arc n pulses a stimulated echo will not 
be create d 43 ' 61 . Nevertheless if a spatial distribution of 
flip angles is present in the sample, stimulated echoes 
will be created due to the deviations from tt pulses as 
shown in refi 43 ' 60 . The long tails previously observed 
in inhomogeneous samples are a consequence of a pulse 
angle distribution among the sample. From Eq.([5]l it can 
be deduced that the long tails observed are due to the 
formation of stimulated echoes and its consequent storage 
of local polarization. For example, the magnetization in 
Ceo was stored as polarization for times two orders of 
magnitude longer than the decoherence time (T 2 he)^ _ — • 

V. CONTROLLED GENERATION OF 
PSEUDOCOHERENCES 

In this section, we introduce a method to artificially 
create, even in homogeneous samples, the conditions un- 
der which the stimulated echoes are formed. The purpose 
of this method is to study in detail the amplitude of the 
polarization stored and the time the polarization can be 
stored as a function of the inhomogeneity added to the 
sample. Sequences can be developed in order to take 
advantage of this stored polarization—. 

The system we have previously studied, Cqo, presents 
an inhomogencously broadened line which satisfies 
T2HE ^ 7T 2 *- Similar situations were observed in all the 
systems that manifest long decay times when CPMG-likc 
sequences are applie d 12 ' 43 ' 44 . 

Here, we use a homogeneous sample to mimic the be- 
haviour observed in inhomogeneous samples by adding 
magnetic field gradients. When a field gradient, G = 
dB z /dz, is applied in the direction of the static mag- 
netic field Boz, the z component of the magnetic field is 
written as: 



B z = B + zG. (9) 

As a consequence of the field gradient presence, the 
spin frequency is modified as follows: 

lo(z) = ujq + ujq(z) (10) 

where log{ z ) = lGz, with 7 the gyromagnctic factor. 

Similarly to Eq.Q, the field gradient produces a 
Hamiltonian linear in spin operators. If a strong enough 
gradient is applied, the main interaction for the spin evo- 
lution is: 

U G = ~fhGY,Itzi- (11) 

The presence of this Hamiltonian term causes an in- 
homogencously broadened line whose width depends on 
the applied field gradient. 

In the present work the sample used as example is poly- 
dimcthylsiloxanc (PDMS)SiS. The characteristic time 
decays for the PDMS arc T 2 * « 0.6 ms and T 2H e ~ 1.8 ms 
indicating a small inhomogeneity in the system. The con- 
dition T2HE ~ 7 Tj found in Ceo is easily satisfied with a 
field gradient of strength around 5 G/cm. All the exper- 
iments were performed at 7 Tesla in a Bruker Avancell. 
A birdcage coil of 10 mm ID and Bruker microimaging 
coils were used. The pulse durations were 12.5 [is for the 
^ and 25.5 fis for the ir pulses. 

To validate the proposed method we applied a (^) — 
t — (71") — t\ — (71") sequence with r = 0, 5 ms and 
t\ = 8 ms with and without field gradient, see Figj6] 
The temporal windows are chosen such that the refocal- 
ization of Hahn and stimulated echoes are produced at 
different times, remember that for perfect tt pulses no 
stimulated echoes are expected. However, even in ab- 
sence of the external field gradient a small stimulated 
echo is observed due to the sample inhomogeneity men- 
tioned before. Notice that its amplitude is approximately 
eight times smaller than the Hahn echo amplitude after 
approximately 4 x T 2 he- Upon application of the external 
field gradient G, the Hahn echo signal is spread into dif- 
ferent frequency components keeping the amplitude con- 
stant, while the stimulated echo increases depending on 
G. It is clear that under the interpulse timings required 
by the CPMG (h=2r) the STE overlaps the Hahn echo. 
By applying proper phase cycling schemes, the discrimi- 
nation of these two echoes becomes possible 4 ^. 

In Fig. [7] the STE is normalized with the amplitude 
of the HE to give its relative contribution as a function 
of the applied gradient field G. Two datasets were ac- 
quired with the proper phase cyclings in order to obtain 
the stimulated and the Hahn echoes independently. It is 
noticeable that the stimulated echo grows quadratically 
with the magnitude of the applied field gradient. 
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FIG. 6: Stimulated echo formation in presence of field gradi- 
ents for the three-pulse sequence shown at the bottom of the 
figure. Notice that while the Hahn echo amplitude remains 
constant the stimulated echo increases with the magnitude of 
the applied field gradient. 



FIG. 8: Magnetization as a function of time for r = 200 fis 
and several field gradient strengths in the CP2 sequence. The 
Hahn echo and CP2 decays are very similar in the absence 
of field gradients. The inset shows the evolution of the mag- 
netization for different values of r when G = 0. Notice the 
change of the trend in the amplitudes between r = 300 and 
r = 400 [is 
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FIG. 7: Stimulated echo amplitude normalized with the cor- 
responding Hahn echo amplitude as a function of the applied 
field gradient. By fitting the data with the function aG + b 
the parameters a = 0.00137 (^) and b — 0.02 are ob- 
tained, straight line in the plot. 



VI. LONG TAILS GENERATED USING 
MAGNETIC FIELD GRADIENTS 

As reviewed in Section III, and Figs. 1-4, we analyzed^, 
two different sequences derived from Carr-Purcell^ and 
Carr-Purcell-Meiboom-Gill^i which produce long decay 
times compared to the Hahn echo in inhomogeneous sam- 
ples, namely CP2 and CPMG1 (see FigH}. In this sec- 
tion, the evolution of the magnetization during these se- 
quences is studied upon variations of external field gra- 
dients and interpulse spacing. In absence of an external 
field gradient the decays observed with the multipulsc 



CPMG1 and CP2 and the Hahn echo sequences in PDMS 
are almost the same. There is however a low contribution 
(c.a. 3%) with decay time of 40 ms that can be observed 
by using a logarithmic scale in the inset of Fig. [H This is 
a consequence of the very small contribution of the stim- 
ulated echo generated by the small inhomogcncity of the 
line even when G = 0. Notice that the amplitudes of the 
tails decrease from r = 100 [is to r = 300 fis, while at 
t = 400 /is this trend reverses. 

In Fig. [5] we show the dependence of the magnetization 
as a function of time for different field gradicntes for the 
CP2 sequence at a fixed value r = 200 /is. As stronger 
field gradients are applied during the complete pulse train 
for both CP2 and CPMG1, bigger amplitudes of the slow 
decaying contributions appear manifesting the storage of 
coherences. 

As observed in C60, the magnetization decay in PDMS 
under gradients manifests two different contributions. 
For short times, the Hahn echo is dominant while for long 
times only the stimulated echo remains, storing polariza- 
tion which is then recovered as transverse magnetization. 
In order to characterize the decays the data were fitted 
with a double exponential, where the short characteristic 
time was fixed to the T2HE decay time. 

M(2t) = A s exp(2r/T 2ffB ) + A t exp(2r/t,) 

The parameters A s , A; and t;, where the subindexes s 
and I stand for short and long decays respectively, were 
calculated by fitting each dataset individually . The pa- 
rameters A; and t; are signatures of the stimulated echoes 
manifested through the long tails. A more detailed anal- 
ysis of the long tails was achieved by discarding the short 
time points and fitting to a single exponential. In Figs.JS] 
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and [TO] the percentage of magnetization stored as polar- 
ization (Ai/A s x 100) and the corresponding values for 
the inverse of the characteristic long decay times t; re- 
spectively are shown as a function of the inter-pulse time 
scaled by T 2 he in the CPMG2 sequence. From Fig[9j it 
can be seen that for the strongest gradient applied and 
for the shortest r, more than 35% of the initial magne- 
tization was stored as polarization and then refocuscd 
during the multipulse sequence. The characteristic de- 
cay time t; of the long tails, is interpreted as a storage 
time and the scaling with T 2 he emphasizes the fact that 
the magnetization can be preserved for times an order of 
magnitude longer than the single-spin decoherence time. 
Under favorable experimental conditions 2t/T2he can be 
decreased to get a better t;/T2HE ratio. This amount can 
also be improved when shorter delays r are experimen- 
tally affordable. 

A remarkable behaviour of 1/t/ vs. t is plotted FigfTQl 
for different values of the gradient G. A change of regime 
occurs at approximately r c = 450 /is corresponding to 
t; f=s 40 ms, the same relaxation time that is obtained in 
the absence of gradients. For values of r shorter than 
t c , the linewidth Av decreases (i.e. the characteristic 
long time tj increases) as the magnitude of the gradient 
increases (see Fig. [TTj) . For times longer than r c , the 
behavior reverses, i.e. t/ decreases as G increases. This 
is more clearly appreciated in Fig llll where the values 
of t; for different values of r are plotted vs. the mag- 
nitude of the gradient G. It is seen that for r < 400/is 
t; increases with the gradient while for t > 500yus t/ de- 
creases with G. This t c (interpulse time 900 /xs) should 
be related to the characteristic fluctuation time of the 
bath t;,, or its spectral density, essentially dominated by 
dipolar interactions partially averaged by motion. In a 
simple picture, one can interpret that a fast interpulse 
rate (1 / (2t c )^>1 / Vfe) produces a decoupling of the sys- 
tem from the bath, leading to narrower lines (longer t; 
values) while the opposite condition leads to a broaden- 
ing of the line. In other terms, one can assume that at 
short values of r (t <C where d stands for the mag- 
nitude of the dipolar interaction, there is an "effective 
dipolar interaction" given by d e ff = d 2 r as the pulses 
act as a breaking time or 1/t as an exchange rate in 
an exchange narrowing process^. In this regime, dipo- 
lar spin dynamics are very slow. When t > 1/d, the 
spin diffusion, which is able to average the deviation of 
the 7r pulses which lead to the formation of stimulated 
echoes, is also effective to blur out the gradients in the 
sample. Then, even with larger gradients, the values of 
tj decrease spoiling the preservation of coherences as po- 
larization. While a deeper analysis of how the dynamical 
variables d, r, and 7G interact is beyond the scope of 
this article, one conclusion is clear: in order to save co- 
herences as polarization the shortest t and the biggest 
G should be used. This can even allow a manipulation 
of nuclear spins with better spatial resolution. As it is 
clear from Eq. ^ , and schematized in FigJSJ the STE se- 
quence is not able to recover the whole coherence, but 
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FIG. 9: Percentage of magnetization stored as polarization 
as function of the inter- pulse time 2r scaled by T2HE in the 
multipulse sequence CP2. Better performance is achieved by 
increasing the field gradients and/or decreasing r. 
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FIG. 10: Characteristic linewidth 1/t; as function of r in the 
multipulse sequence CP2. Note the crossing point at approx- 
imately 450 /is 



at most 1/2 of its initial amplitude. Consequently, this 
strategy is not useful to make "quantum memories" (fi- 
delity below 0.66). Nevertheless, it is potentially useful 
to make read out memories. Moreover, in a different con- 
text, Kobayashi et al<^ have designed variations of the 
STE sequence that overcome the 1/2 factor and recover 
a much bigger fraction of the coherences (or transverse 
magnetization) . 

Finally, we address the possibility of performing exper- 
iments with multipulse sequences that do not produce 
coherence storage as polarization. A four phase cycle 
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FIG. 11: Characteristic long decay rate as a function of the 
field gradient. Different values of r, which scale the dipolar 
dynamics (d c g, ~ d 2 r), from top to bottom, 200, 300, 400, 
500, 600, 700, 800, 900 and 1000 fj,s are displayed. Note the 
change of slope between 400 and 500 /is. This is emphasized in 
the inset where the slopes are plotted as a function of r. The 
change of sign identifies the critical value. The decay of ti = 
40ms for G — is due to line inhomegeneities and shimming 
imperfections 



sequence that resembles the MLEV-4^: 

CPMG4 : (-) — {t — iTy — t — echo 
V 2 / x 



T 
T 
T 



t — echo 

— t — echo 

— t — echo ) r 



(12) 



has been empirically shown to be the most suitable 
for preventing unexpected artifacts in a T2 measure- 
ment^. 

In Fig. [T2]the results of a Hahn echo and CPMG4 mea- 
surements with and without field gradient are plotted. 
It can be observed that the magnetization decay agrees 
perfectly well for the three cases. The reason why this 
sequence is the appropriate to measure T2 can be un- 
derstood from the stimulated echo perspective presented 
previously^. The CPMG4 sequence does not accumu- 
late magnetization contributions from the constructive 
interference between normal and stimulated echoes, and 
this is clearly manifested in the results shown in Fig. 1121 
Even in the presence of an inhomogeneous field, long tails 
are not observed when this sequence is applied. 



VII. CONCLUSIONS 

In this article, we have reviewed the conditions for the 
formation of stimulated echoes and the consequent gen- 
eration of long tails in the magnetization decay. Besides, 
we report an original set of experiments using field gradi- 
ents in homogeneous samples that extend and validates 



FIG. 12: Magnetization observed by applying the CPMG4 
with gradient off and on compared to the Hahn Echo decay. 



the previous findings and their interpretation. Moreover, 
they suggest a controlled way to preserve coherences as 
polarization that may account for more than 1/3 of the 
original signal. The long tails observed when multipulse 
sequences as CPMG are applied to inhomogeneous sam- 
ples were assigned to polarization that encodes quantum 
coherences. In this way the coherences are stored for 
times much longer than T2he- It has been shown that 
the responsible of the long tails are the stimulated echoes, 
formed due to the inhomogeneity of the sample. The re- 
quirements to observe these slow magnetization decays 
are (1) an rf field inhomogeneity or a highly inhomo- 
geneous line able to produce different tilting angles in 
different sites of the sample and (2) non-effective flip- 
flop interactions. Under these conditions the differences 
in resonance frequencies are larger than the dipolar cou- 
plings, making ineffective the flip-flop mechanism. Con- 
sequently, the different deviations from 7r pulses (regard- 
less of its origin) cannot be averaged out by dipolar inter- 
actions. For 2r « T2HE the amplitude of the coherence 
stored is less than 5%. In this condition, the time be- 
tween pulses is enough for the dipolar interaction to be 
operative. 

In this paper we showed that the stimulated echoes pre- 
serve coherent states frozen in the direction of the static 
magnetic field, then bring them back to the plane where 
they are refocused and observed. It is important to re- 
mark that the magnetization decay is very long because 
the coherences are stored as polarization. So, these long 
decays should not be interpreted as a long decoherence 
time during which quantum operations can be performed. 
However, this strategy can be exploited to design mem- 
ories in nuclear or electron spin systems. Then, upon re- 
covery of the stored coherences, interactions among the 
whole system can be turned on. The concept introduced 
in this paper was indeed indirectly exploited in Ref^S. 
There, coherences induced in electron spins were trans- 
ferred to a coherent state between nuclear spins, stored 
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by applying the CPMG1 sequence and transferred back 
to electron spins. 

The method that mimic the conditions to induce the 
stimulated echoes for homogeneous samples in a con- 
trolled way is based on the application of field gradients 
that produce a frequency distribution across the sample. 
This, in turn, generates a pulse angle distribution that 
sets the conditions under which the stimulated echoes ap- 
pear. An important fact deduced from these experiments 
is that the origin of the inhomogeneity of the samples is 
not relevant for the observation of long tails. Whether 
the inhomogeneity arises from a random distribution of 
nuclear spins (as occurs in 29 Si) or on a linear dependent 
gradient, like the one used in the present experiments, the 
consequences are the same. Moreover, the use of suitable 
field gradients could allow a more selective spatial ma- 
nipulation of the memories. This method was applied to 
a PDMS sample, for which we were able to perform a 
complete quantification of the procedure, i.e. the stor- 
age time and ratio of coherence survival. There, it was 
possible to detect a transition where the effective dipolar 
interaction d 2 r becomes of the order of the main frequen- 
cies in the bath's spectral density. This imposes a limit in 



the interpulsc spacing in which coherences can be stored. 
Under the condition r < l/d higher polarization storage 
and longer decay times are achieved by using stronger 
gradients. 

Finally, we tested the CPMG4 sequence under field 
gradients that reproduce inhomogeneous properties. The 
mechanism behind the correct operation of this pulse se- 
quence can be understood under the stimulated echoes 
interferences. The true single-spin T2 decay time can be 
measured under this pulse sequence either for homoge- 
neous or inhomogeneous samples. Thus, we verified the 
empirically known robustness of this sequence which, re- 
gardless of the gradient, yields the same results as the 
Hahn echo. 
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